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CALIBRATED  L-BAND  TERRAIN  MEASUREMENTS 
AND  ANALYSIS  PROGRAM  -  RESULTS 

1 

INTRODUCTION 

This  report  presents  the  results  obtained  during  the  three  years  of 
the  Calibrated  L-band  Terrain  Measurements  and  Analysis  Program  con¬ 
ducted  by  the  Environmental  Research  Institute  of  Michigan  (ERIM)  for 
the  Rome  Air  Development  Center  of  the  USAF.  The  overall  goal  of  this 
program  was  to  provide  L-band  radar  clutter  data  from  a  variety  of 
terrain  types  and  sea  states  and  to  begin  to  evaluate  models  which 
describe  the  radar  clutter.  The  approach  adopted  for  this  program  was 
to  utilize  airborne  synthetic  aperture  radar  (SAR)  (Rawson,  et  al., 
1985)  Imagery  as  a  source  for  the  L-band  clutter  data.  The  SAR  Imag¬ 
ery  used  during  this  program  was  provided  courtesy  of  several  differ¬ 
ent  U.S.  Government  sponsors.  The  primary  data  sets  were  collected 
over  central  North  Carolina  during  a  U.S.  Geological  Survey  Radar  data 
collection  (Kaslschke,  1985;  Kover  and  Jones,  1985)  and  over  Long 
Island,  New  York  during  the  ONR-sponsored  SARSEX  experiment 
(Kaslschke,  et  al.,  1985).  Additional  SAR  data  was  collected  over  the 
Marginal  Ice  Zone  during  the  ONR-sponsored  MIZEX  '84  experiment  (MIZEX 
Group,  1986). 

This  program  was  conducted  In  two  distinct  phases.  During  Phase 
I,  four  distinct  tasks  were  performed:  (1)  collection  of  SAR  data  and 
associated  ground  truth  Information;  (2)  processing  and  calibration  of 
SAR  data  of  selected  test  areas;  (3)  reduction  of  ground-truth  data; 
and  (4)  selection  of  radar  surface  scattering  models.  The  Phase  I 
activities  are  summarized  by  Larson,  et  al.  (1986).  During  Phase  II, 
the  following  tasks  were  performed:  (1)  extraction  of  calibrated  L- 
band  radar  cross  section  measurements  from  the  SAR  data; 
(2)  evaluation  and  modeling  of  the  Image  intensity  distributions 
observed  on  SAR  Imagery;  and  (3)  comparison  of  the  radar  scattering 
coefficients  obtained  In  Phase  1  to  the  ground-truth  data  via 


deterministic  models  and  statistical  analyses.  The  procedures  used 
to  analyze  the  SAR  Imagery,  along  with  the  results  of  these  analyses 
are  presented  In  this  report. 

Including  this  Introduction,  this  report  Includes  five  chapters. 
Chapter  2  presents  a  summary  of  the  radar  scattering  coefficients 
extracted  during  this  program.  In  addition  to  the  L-band  SAR  data,  a 
limited  amount  of  X-band  and  C-band  radar  scattering  coefficient 
measurements  were  made  available  to  this  program.  These  measurements 
are  summarized  In  this  chapter.  Chapter  3  of  this  report  discusses 
the  modeling  of  the  Intensity  distributions  observed  on  SAR  data. 
Chapter  4  presents  the  results  of  the  comparison  of  the  SAR  radar 
cross  sections  to  the  ground-truth  data.  Finally,  Chapter  5 
discusses  the  conclusions  and  recommendations  from  this  program. 
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2 

L-8ANB  RADAR  SCATTERING  COEFFICIENTS 

In  this  chapter,  the  radar  scattering  coefficient  (a0)  measure¬ 
ments  derived  from  the  various  SAR  data  sets  are  presented.  Calibra¬ 
tion  of  the  SAR  system  Is  described  In  the  Interim  report  (Larson,  et 
al  1986).  Because  the  size  of  the  individual  test  sites  within  the 
several  different  data  sets  varied  considerably,  different  data  ex¬ 
traction  and  display  techniques  were  employed.  These  techniques  will 
be  discussed  In  the  different  sections  of  this  chapter. 

In  this  chapter,  we  will  present  a0  measurements  from  four  dif¬ 
ferent  test  areas: 

1  Long  Island,  New  York  (SARSEX) 

2.  Duke  Forest,  North  Carolina  (USGS) 

3.  Atlantic  Ocean  (SARSEX),  and 

4.  Fram  Strait,  Greenland  Sea  (MIZEX). 

2.1  LONG  ISLAND,  NEW  YORK  TEST  SITES 

The  Long  Island,  New  York  test  sites  were  Imaged  during  a  number 
of  different  passes  during  the  SARSEX  experiment.  The  test  sites  used 
for  this  study  were  located  In  and  around  the  Peconlc  River  Airport, 
which  Is  owned  and  operated  by  Grumman  Aerospace  Corporation.  For 
this  analysis,  we  will  present  data  collected  during  six  different  SAR 
passes  on  two  separate  dates  at  different  Incidence  angles  and  polari¬ 
zations.  Specifically,  we  will  present  data  from  four  passes  where 
HH-polarlzatlon  data  were  collected  and  data  from  two  passes  where  VV- 
polarlzatlon  data  were  collected.  Figures  1  and  2  present  examples  of 
the  L-band  HH-  and  VV-polarlzed  SAR  Imagery  collected  over  the  Long 
Island  test  sites. 

Radar  scattering  coefficient  (a0)  measurements  were  extracted  from 
22  different  test  sites.  These  test  sites  are  summarized  In  Table  1 


3 


Figure  l.  L-band  (HH)  SAR  Imagery  of  the  Long  Island  Test 
Sites  Collected  During  SARSEX-10,  Pass  1,  7 
October  1984 
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Figure  2.  L-band  (W)  SAR  Imagery  of  the  Long  Island  Test 
Sites  Collected  During  SARSEX-10,  Pass  3,  7 
October  1984 
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TABLE  1 


Site 


HUB m 


LONG  ISLAND  TEST  SITES 


Designation 


Site  Description 


A 

B 

C 

D 

E 

0 

8 

3 

F 

H 

I 

J 

K 

RF1 

RF2 

RF3 

L 

M 

N 

0 

P 

Q 

R 

W 

z 


Strawberry  Field 
Orchard 

Smooth,  Bare  Field 

Smooth  Field  -  Near  South  Rd. 

Potatoes  -  South  Side 

Potatoes 

Potatoes 

Potatoes 

Potatoes 

Grass 

Corn  Stalks 

Grass 

Cabbage 

Grass  Field  (VOR,  Reflectors) 
Grass  Field 
Grass  Field 
Grass  Field 

Short  Pike,  Trees  -  Weeds 
Smooth,  Bare  Field 
Smooth,  Bare  Field 
Smooth,  Bare  Field 
Smooth,  Bare  Field 
Sod 
Sod 
Sod 

Bare,  Smooth  Field,  N-S  Furrows 
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and  their  locations  presented  In  Figure  3.  Figure  4  presents  selected 
surface  photographs  Illustrating  the  different  types  of  ground  cover 
present  at  the  Long  Island  test  sites.  (Other  photographs  included  In 
Larson,  et  al  1986). 

The  a0  measurements  were  generated  by  averaging  a  40  by  40  pixel 
subset  from  the  Individual  test  sites.  The  ERIM  SAR  Image  processor 
generates  pixels  with  a  dimension  of  1.5  by  1.44  m.  The  3  m  resolu¬ 
tion  of  the  SAR  data  Is  therefore  oversampled  by  a  factor  of  2  In  both 
the  range  and  azimuth  dimensions.  In  order  to  estimate  the  number  of 
Independent  looks,  N,  for  the  o°  measurements,  we  must  divide  the 
number  of  pixels  averaged  by  2  In  both  the  range  and  azimuth  di¬ 
mensions  (Kaslschke,  et  al.,  1987).  Each  a0  value  thus  was  generated 
using  400  Independent  samples  (20  by  20  pixels),  resulting  In  a  90% 
confidence  Interval  due  to  speckle  or  fading  of  +  0.35  dB,  assuming 
that  radar  speckle  has  a  chi-square  distribution  (Ulaby,  et  al., 
1982). 

Another  source  of  uncertainty  In  the  <7°  measurements  Is  the  within 
field  variation  present  In  the  data.  To  measure  this  uncertainty,  six 
different  a0  measurements  were  obtained  from  three  of  the  test  sites 
(F,  G  and  H)  at  four  separate  Incidence  angles  using  the  L-band  (HH) 
data  sets.  These  measurements  are  summarized  In  Table  2.  The  maximum 
standard  deviation  obtained  was  0.5  dB.  We  can  define  a  standard 
error  of  the  mean  (SE)  as 


SE  =  SD/Nl/2  (!) 

where  SO  Is  the  standard  deviation  and  N  Is  the  number  of  samples 
averaged.  Doing  so  results  In  a  maximum  standard  error  of  +  0.3  dB. 
Combining  the  uncertainty  due  to  radar  fading,  with  the  maximum  uncer¬ 
tainty  due  to  the  within  field  variability  results  In  a  maximum  +  0.5 
dB  error  bound  for  the  o°  measurements.  In  Figure  5,  we  plot 
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Figure  5.  Plot  of  Average  L-band  (HH)  o°  Values  for  Three 
Test  Sites 
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the  cr*  measurements  presented  In  Table  2.  We  can  see  that  for  all  3 
test  sites,  there  Is  a  clear  decrease  In  a*  as  Incidence  angle 
Increases.  In  these  examples,  the  a*  values  for  the  standing  corn 
stalks  are  clearly  5  dB  higher  than  the  other  two  test  sites,  whose 
a *  values  are  essentially  equivalent. 

The  a*  values  for  the  14  test  sites  are  summarized  In  Figure  6. 
From  these  graphs,  we  can  note  the  following  trends.  In  general,  for 
L-band  (HH)  radar  data,  as  Incidence  angle  Increases,  a*  decreases. 
This  trend  was  observed  for  all  data  sets.  The  trend  for  the  L-band 
(VV)  data  was  not  as  clear,  with  c*  both  increasing  and  decreasing  as 
Incidence  angle  Increases. 

Figure  /  presents  a  composite  plot  of  the  various  fields  compared 
to  one  another.  In  this  plot,  we  combined  all  similar  test  sites  to 
obtain  an  average  a*  signature  for  a  particular  terrain  type.  In 
this  plot,  we  compare  only  the  L-band  (HH)  data  sets.  The  highest  a* 
values  were  obtained  for  the  orchard  and  corn  stalks,  while  the 
lowest  were  obtained  for  the  sod  and  potato  fields.  There  Is  clearly 
a  10  dB  difference  between  the  darkest  and  brightest  fields  at  all 
Incidence  angles. 

2.2  DUKE  FOREST  TEST  SITES 

In  early  April  1984,  during  a  data  collection  program  conducted  for 
the  U.S.  Geological  Survey  (Kover  and  Jones,  1985),  multifrequency, 
multipolarization  SAR  data  were  collected  over  the  Duke  Forest  using 
the  ERIM/CCRS  CV580  SAR  System  (see  Figures  8  and  9).  In  this  study 
we  utilize  L-band  (VV)  data  collected  on  8  April  1984  and  X  and  C- 
band  (VV)  data  collected  on  10  April  1984.  The  deciduous  trees  In 
the  study  area  were  at  a  stage  where  most  of  the  buds  had  burst,  but 
no  significant  leaf  flushing  or  new  stem  growth  had  occurred.  The 
weather  conditions  were  Identical  on  each  day  of  the  SAR  data 
collection,  with  overcast  skies  and  periods  of  light  rain. 
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The  tree  stands  In  the  Duke  Forest  range  In  size  from  1  to  >  50 
hectares,  and  have  been  In  a  "forest  condition"  for  more  than  50 
years.  What  makes  this  a  particularly  valuable  test  area  for  the 
evaluation  of  SAR  Imagery  for  forestry  studies  Is  that  research  rec¬ 
ords  for  the  stands  within  the  forest  are  extensive  (Edeburn,  1981), 
with  historical  Information  on  tree  growth,  stand  age,  stand  density, 
stand  mortality,  etc.,  In  existence  for  hundreds  of  permanent  plots 
throughout  the  forest.  In  addition,  the  Duke  Forest  Is  the  site  of 
numerous  ongoing  and  proposed  forestry  and  ecological  studies  (see, 
e.g.,  Christensen  and  Peet,  1984). 

For  the  present  analysis,  a  2.9  by  6.1  km  area  surrounding  the 
Blackwood  Division  of  the  Ouke  Forest  was  utilized  as  the  primary  test 
area.  Figure  10  Is  a  mosaic  of  1:9600  scale  black  and  white  aerial 
photographs  which  contains  the  Blackwood  Division.  The  aerial  photo¬ 
graphy  was  collected  In  January,  1985,  nine  months  after  the  SAR  data 
collection.  Figure  11  presents  a  vegetation  map  of  the  Blackwood 
Division  of  the  Duke  Forest  which  was  generated  using  a  combination  of 
existing  forest  cover  maps  (Edeburn,  1981),  ground  truth  collected  at 
the  time  of  the  SAR  overflights,  and  the  aerial  photography  presented 
in  Figure  10.  A  point  of  reference  for  all  the  Images  and  maps  pre¬ 
sented  In  this  paper  Is  the  large  grass  field  (site  B9)  located  In  the 
center  of  a  recent  clearcut  area  of  the  Blackwood  Division. 

A  secondary  test  area  outside  of  the  Blackwood  Division,  which 
contained  a  forest  stand  located  along  a  stream  which  was  flooded 
during  the  SAR  overflights,  was  also  used  in  this  study.  An  enlarge¬ 
ment  of  the  SAR  Imagery  from  this  secondary  test  area  Is  presented  in 
Figure  12. 

Twenty-two  test  sites  were  Identified  within  the  two  test  areas 
(see  Table  3  and  Figure  10).  These  test  sites  can  be  divided  into  the 
following  basic  categories: 

Young  Pine  (Figure  13a):  A  three  year  old  stand  of  loblolly  pine 

(P.taeda),  1.2  m  In  height  (site  B8). 
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Figure  10.  Aerial  Photographic-Mosaic  Collected  over  the 
Blackwood  Division  of  the  Duke  Forest 
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C-Band 


L-Band 

Figure  12.  Digitally-Processed  X-,  c-  and  L-band  SAR  Data 
(W  polarization)  Collected  over  a  Flooded 
Forest  Stand  North  of  the  Blackwood  Division  of 
the  Duke  Forest 
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TABLE  3 

DUKE  FOREST  TEST  SITES 


CONFIGEROUS  SPECIES  DECIDOUS  SPECIES 

BASAL 


TREE 

STEMS/ 

DBH 

AREA 

HEIGHT 

STEMS/ 

DBH 

BASAL 

HEIGHT 

SITE 

TYPE 

AGE 

ha 

(cm) 

(sq  m/ha) 

(m) 

ha 

(cm) 

AREA 

(■) 

B1 

0 

U 

7 

36 

1 

21 

640 

17 

27 

22 

B2 

SO 

U 

99 

30 

8 

26 

388 

19 

18 

B3 

C 

U 

121 

37 

18 

27 

B4 

S 

60 

121 

34 

11 

31 

1247 

11 

18 

B5 

AO 

U 

10 

54 

2 

31 

472 

19 

22 

B6 

LE 

U 

106 

36 

12 

25 

810 

19 

22 

B7 

C 

U 

10 

41 

1 

27 

699 

18 

30 

B8 

L 

3 

8 

1 

69 

G 

- 

BIO 

L 

50 

247 

35 

32 

26 

378 

15 

9 

BU 

C 

U 

10 

41 

1 

24 

566 

28 

27 

24 

B12 

D 

U 

104 

27 

7 

23 

299 

31 

30 

24 

B13 

D 

U 

20 

B14 

L 

60 

170 

39 

21 

30 

729 

16 

24 

B15 

L 

40 

454 

31 

37 

27 

17 

41 

3 

B16 

L 

40 

437 

28 

28 

24 

284 

10 

3 

B18 

L 

30 

410 

22 

18 

20 

B19 

L 

30 

963 

24 

46 

B20 

L 

30 

531 

25 

28 

22 

79 

21 

5 

B21 

E 

U 

442 

22 

27 

28 

B22 

L 

40 

358 

31 

28 

27 

823 

L 

30 

390 

25 

21 

22 

72 

20 

3 

824  SSF  U 
625  SSD  U 

TREE  TYPE  KEY 

AO  -  SWEETGUM/YELLOW  POPLAR 

C  -  WHITE  OAK/REO  OAK/BLCK  OAK 

D  -  BLACKJACK  OAK/POST  OAK 

E  -  MIXED  HARDWOOD 

L  -  LOBLOLLY  PINE 

LE  -  LOBLOLLY  PINE/MIXED  HARDWOOD 

S  -  SHORTLEAF  PINE 

SO  -  SHORTLEAF  PINE/OAK 

SSD  -  SYCAMORE/SWEETGUN/LOBLOLLY-DRY  SITE 

SSF  -  SYCAMORE/SWEETGUM/LOBLOLLY-FLOODED  SITE 
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Nature  Pine  (Figure  13b):  Thirty  to  forty  year  old  loblolly  pine 
stands,  with  virtually  no  other  tree  species  In  the  overstory  or 
understory  (site  823). 

Old  Pine  (Figure  13c):  50  to  60  year  old  stands  of  loblolly  pine 
and  shortleaf  pine  (P.  echlnata),  where  mortality  In  the  overstory 
has  allowed  the  Invasion  of  deciduous  species  (site  B14). 

Nixed  Pine/Hardwood  (Figure  13d):  Stands  where  both  pine  and 
hardwoods  are  dominant  overstory  species  (site  B6). 

Hardwood  Stands  (Figure  13e):  Stands  where  hardwood  species  domi¬ 
nate  the  overstory.  On  drier  sites,  pure  and  mixed  stands  of  oak 
prevail  (Q.  alba,  Q.  rubra,  Q.  velutlna,  Q.  marllandlca,  Q.  prlnus 
and  Q.  stellata),  while  on  molster  sites,  yellow  poplar  (L.  tullp- 
Ifera)  and  sweetgum  (L.  styraclflua)  are  the  dominant  species 
(site  B7) . 

Flooded  Stand  (Figure  13f):  Spring  rains  resulted  In  one  stand  of 
sycamore  (P.  occidentalism  and  sweetgum  being  flooded  at  the  time 
of  the  SAR  data  collection.  Although  the  picture  of  this  area  in 
Figure  13f  was  collected  at  a  time  when  the  site  was  not  In  a 
flooded  state,  the  water  marks  on  the  tree  trunks  clearly  Indicate 
that  this  stand  Is  flooded  at  some  stage  during  the  year  (site 
B24) . 

Grass  Field  (Figure  13g):  This  test  site  containing  short  (5  cm 
In  height)  grass  was  included  as  a  reference  for  the  other  test 
sites  (site  B9). 

The  research  records  from  the  Duke  Forest  were  reviewed,  and  a  set 
of  tree  parameters  obtained  for  each  test  site.  These  parameters 
Include: 

Age:  The  average  age  (usually  rounded  to  the  nearest  ten  years) 
of  the  trees  within  the  stand. 

Items/hectare:  The  average  number  of  trees  per  hectare  greater 
than  10  cm  In  diameter  within  the  stand. 

Diameter  at  breast  height  (DBH):  The  average  diameter  of  all 
trees  within  the  stand  measured  at  a  height  of  1.3  nr  above  the 
ground  level. 

Basal  area:  The  total  area  (In  m2)  of  all  tree  trunks  (at  1.3  m 
above  the  ground)  within  the  forest  stand. 

Height  (m):  The  average  height  of  the  trees  within  the  stand. 


38 


These  parameters  are  summarized  for  the  22  test  sites  In  Table  3. 
Note  that  the  data  presented  In  Table  3  has  been  divided  Into  two 
groups  (coniferous  and  deciduous  species)  for  each  test  site. 

Estimates  of  the  radar  scattering  coefficient  (*°)  were  generated 
for  the  22  test  sites  using  the  average  value  from  a  40  by  40  pixel 
subset  extracted  from  the  calibrated  SAR  Images. 

Table  4  summarizes  the  a0  values  for  the  22  test  sites.  The  Inci¬ 
dence  angle  ($j)  for  the  test  sites  Is  57°,  except  for  sites  B24  and 
B25,  where  It  Is  46°.  In  Table  4  the  test  sites  have  been  grouped  by 
tree  type  and  age.  For  analysis  purposes,  the  22  test  sites  were 
organized  Into  13  categories.  The  average  a0  values  for  these  13 
categories  are  summarized  in  Table  5  and  presented  In  Figure  14.  The 
upper  and  lower  error  bounds  are  also  presented  In  Figures  14a-14c. 
We  define  the  upper  bound  (UB)  being  equal  to  the  lower  bound  (LB)  as 

UB  *  LB  -  [r2  +  Sf2  +  SE2] 1/2  (2) 

where  r  is  the  relative  calibration  of  the  SAR  (1  dB),  Sf  Is  the  var¬ 
iability  due  to  speckle  (.35  dB)  and  SE  Is  the  standard  error  of  the 
mean  defined  In  Eq.  (1).  Figure  14d  presents  the  average  a0  values 
for  all  three  radar  frequencies  plotted  on  the  same  graph  for  compari¬ 
son  purposes. 

From  Tables  4  and  5  and  Figure  14,  we  can  make  several  observa¬ 
tions.  We  first  note  that  while  there  Is  a  significant  difference 
between  the  dry  and  flooded  sycamore  stands  In  the  L-band  (3.5  dB)  and 
C-band  (3.0  dB)  Imagery,  there  Is  considerably  less  difference  In  the 
X-band  data  (0.9  dB) .  The  observation  at  L-band  Is  consistent  with 
studies  conducted  with  satellite  SARs  (Krohn,  et  al . ,  1983;  Ormsby,  et 
al . ,  1985;  Lyon  and  McCarthy,  1981).  At  C-band,  no  previous  SAR  data 
exists  where  flooded  forests  had  been  Imaged.  However,  previous  air¬ 
craft  X-band  SLAR  Imagery  collected  over  a  flooded  forest  canopy 


TABLE  4 

RADAR  SCATTERING  COEFFICIENTS  FOR  DUKE  FOREST  TEST  SITES 


SCATTERING  COEFFICIENTS 

TREE 


SITE 

TYPE 

AGE 

X-BAND 

C-BAND 

L-BAND 

B9 

G 

» 

-5.14 

-5.50 

-16.34 

68 

L 

3 

-1.33 

-3.30 

-13.39 

B19 

L 

30 

-3.11 

-3.97 

-10.98 

B20 

L 

30 

-1.75 

-5.17 

-11.61 

B23 

L 

30 

-1.26 

-  .41 

-11.58 

B15 

L 

40 

-2.88 

-3.51 

-12.12 

B16 

L 

40 

-3.00 

-3.58 

-10.80 

B22 

L 

40 

-4.43 

-3.05 

-11.02 

BIO 

L 

50 

-4.43 

-3.05 

-11.02 

B4 

S 

60 

-2.22 

-4.51 

-13.05 

814 

L 

60 

-1.84 

-4.76 

-11.58 

B2 

SO 

U 

-  .64 

-3.26 

-11.36 

B6 

LE 

U 

-  .74 

-3.67 

-11.29 

B1 

D 

U 

-  .01 

-2.73 

-11.03 

B12 

D 

U 

-3.13 

-3.47 

-11.21 

B13 

D 

U 

-  .73 

-3.23 

-11.67 

B3 

C 

U 

-  .64 

-3.60 

-11.41 

B7 

C 

U 

-  .18 

-2.56 

-10.62 

Bll 

C 

U 

-3.15 

-3.55 

-10.97 

B5 

AO 

U 

-1.30 

-2.77 

-11.28 

B21 

E 

U 

-  .47 

-4.06 

-11.60 

B24 

SSF 

U 

-  .95 

-1.47 

-  8.12 

B25 

SSD 

-1.78 

-5.14 

-11.62 

TREE  TYPE  KEY 

AO  -  SWEETGUM/ YELLOW  POPLAR 

C  -  WHITE  OAK/RED  OAK/BLCK  OAK 

D  -  BLACKJACK  OAK/POST  OAK 

E  -  MIXED  HARDWOOD 

L  -  LOBLOLLY  PINE 

LE  -  LOBLOLLY  PINE/MIXED  HARDWOOD 

S  -  SHORTLEAF  PINE 

SO  -  SHORTLEAF  PINE/OAK 

SSD  -  SYCAMORE/SWEETGUN/LOBLOLLY-DRY  SITE 

SSF  -  SYCAMORE/SWEETGUM/LOBLOLLY-FLOODED  SITE 
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TABLE  5 

MEANS  AND  STANOARO  DEVIATIONS  OF  DUKE  FOREST 
RADAR  CROSS  SECTION  VALUES 


SITE 

N 

X-BAND 

C-BAND 

L-BAND 

GRASS  FIELD 

1 

-9.11 

-0.47 

-16.34 

3  YR  PINE 

1 

-5.31 

-7.28 

-13.39 

30  YR  PINE 

3 

MEAN 

SD 

-6.07 

.95 

-7.21 

2.60 

-11.39 

.29 

40  YR  PINE 

3 

MEAN 

SD 

-6.62 

.66 

-8.10 

.88 

-11.67 

.62 

50  YR  PINE 

1 

-8.36 

-6.98 

-11.02 

60  YR  PINE 

2 

MEAN 

SD 

-6.19 

.23 

-8.79 

.22 

-12.31 

.73 

HARDWOOD/PINE  MIX 

3 

MEAN 

SD 

-4.87 

.06 

-7.55 

.28 

-11.32 

.03 

CHESTNUT  OAK 

3 

MEAN 

SD 

-5.37 

1.30 

-7.22 

.22 

-11.30 

.27 

MIXED  OAK 

3 

MEAN 

SD 

-5.35 

1.29 

-7.22 

.52 

-11.00 

.32 

YELLOW  POPLAR/SWEETGUM 

1 

-5.39 

-6.86 

-11.29 

MIXED  HARDWOODS 

1 

-4.68 

-8.26 

-11.60 

DRY  SYCAMORE 

1 

-6.37 

-9.72 

-11.62 

FLOODED  SYCAMORE 

1 

-5.60 

-6.12 

-  8.12 
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SCATTERING  COEFFICIENT  (<J8)  SCATTERING  COEFFICIENT  (<JB) 


SCATTERING  COEFFICIENT  (<3B)  SCATTERING  COEFFICIENT  (dB) 


(Ormsby,  et  al.,  1985)  revealed  no  significant  difference  from  adja¬ 
cent  unflooded  canopies. 

The  second  observation  we  can  make  from  the  SAR  forestry  data 
concerns  the  ranges  of  radar  scattering  coefficients  observed  over 
different  forest  stands.  In  this  analysis,  we  consider  only  those 
sites  from  the  same  Incidence  angle  range  with  full  canopies  (l.e., 
all  sites  In  Table  4  except  B8,  B9,  B24,  B25) .  If  we  consider  the 
average  values  presented  In  Table  5,  the  range  of  observed  o°'s  at  X- 
band  Is  3.4  dB,  at  C-band,  2.3  dB,  and  L-band,  1.3  dB.  If  we  con¬ 
sider  all  the  values,  as  presented  In  Table  4,  then  the  range  of  scat¬ 
tering  coefficients  for  mature,  dry  forests  Is  4.2  dB  at  X-band, 
4.8  dB  at  C-band  and  2.8  dB  at  L-band.  Thus,  the  highest  range  of 
radar  responses  for  forest  canopies  occurs  at  the  higher  radar  fre¬ 
quencies  (X  and  C-band). 

2.3  ATLANTIC  OCEAN  TEST  SITES 

A  considerable  amount  of  calibrated  L-band  (HH)  SAR  Imagery  was 
collected  over  open  ocean  regions  during  the  1984  SARSEX  experiment. 
In  this  section,  we  will  present  L-band  scattering  coefficients  from 
three  passes  collected  during  this  mission. 

Computer  programs  have  been  developed  to  extract  calibrated  radar 
scattering  cross-section  (a)  or  scattering  coefficient  a0)  measure¬ 
ments  from  airborne  SAR  Imagery.  The  SAR  calibration  algorithm  devel¬ 
oped  by  Larson,  et  al.  (1987)  has  been  coupled  with  an  algorithm  which 
allows  for  scanning  along  any  track  In  any  direction  within  the  SAR 
scene.  This  computer  program  allows  for  the  operator  to  select: 
(1)  scan  direction,  (2)  number  of  pixels  for  averaging  In  the  cross¬ 
track  direction,  and  (3)  number  of  pixels  for  averaging  In  the  along- 
track  direction.  The  algorithm  outputs  a  plot  with  either  slant  or 
ground  range  as  the  X-axis  and  relative  or  absolute  a  or  a0  as  the  Y- 
axls. 


The  tracks  of  the  surface  vessels  collecting  the  ocean  measure¬ 
ments  during  SARSEX  were  always  oriented  perpendicular  to  the  propaga¬ 
tion  direction  of  the  Internal  waves.  Thus,  the  scans  extracted  from 
the  SAR  Images  were  oriented  so  they  were  parallel  to  the  ship's 
track,  with  the  center  of  the  scan  coincident  with  the  ship's  position 
in  the  SAR  Image.  To  generate  these  scans,  a  window  of  150  m  In  the 
cross-track  direction  and  22  m  In  the  along-track  direction  was  used. 
The  number  of  pixels  or  Independent  samples  used  to  generate  each  o° 
measurement  leads  to  a  maximum  +  0.55  dB  uncertainty  (90%  confidence 
Interval)  In  the  resultant  a0  value  due  to  Image  speckle  or  fading 
(Ulaby,  et  al.f  1982). 

The  data  used  In  this  analysis  was  collected  during  three  separate 
passes:  SARSEX-4,  Pass  5;  SARSEX-5,  Pass  9;  and  SARSEX-8,  Pass  4. 
Figure  15  presents  the  a0  scans  for  these  three  passes.  The  Incidence 
angle  ranges  for  these  data  sets  was  from  25°  to  56°.  We  can  see  that 
the  overall  a0  level  for  the  three  passes  was  fairly  constant  with 
respect  to  wind  speed,  but  there  is  a  noticeable  fall-  off  In  o°  as 
Incidence  angle  decreases.  This  fall-off  Is  approximately  0.25  dB  per 
degree  In  Incidence  angle.  The  observed  fall -off  matches  model  pre¬ 
dicted  values  quite  well  (Lyzenga  and  Bennett,  1987;  Kaslschke,  et 
al.,  1987). 

2.4  FRAM  STRAIT  ICE  TEST  SITES 

L-band  (HH)  radar  scattering  coefficient  (o°)  values  were  extract¬ 
ed  from  airborne  SAR  Imagery  collected  during  the  MIZEX  experiment. 
In  order  to  calibrate  this  data  set,  the  L-band  SAR  Image  Intensities 
were  compared  to  helicopter-borne  scatterometer  measurements  collected 
near-col ncldently  with  the  SAR  data  (Onstott  and  Larson,  1986). 

Table  6  summarizes  the  L-band  (HH)  o°  values.  (Burns,  1987)  All 
data  were  collected  at  an  Incidence  angle  of  35°.  Note  that  in  this 
case,  the  o°  value  from  the  open  water  test  site  Is  -32.5  dB, 
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TABLE  6 

SUMMARY  OF  RADAR  SCATTERING  COEFFICIENT  (<7°) 
VALUES  FOR  THE  FRAM  STRAIT  TEST  SITE 


TEST  SITE 
Open  Water 
Brash  Ice 

Thin  First  Year  Ice 
Medium  First  Year  Ice 
Heavy  Snow  Cover 
Pressure  Ridge 
Multiyear  Ice  with  Snow 
Multiyear  Ice  with  Melt 


0  VALUE  (dB) 
-32.5 
-21.5 
-25.0 
-25.0 
-27.5 
-25.5 
-26.0 

Pools  -22.5 
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approximately  15-20  dB  lower  than  the  values  observed  over  the  SARSEX 
test  areas.  Surface  observations  of  the  open  water  areas  In  the  HIZEX 
scenes  showed  that  the  water  surface  was  extremely  calm,  with  no  sur¬ 
face  capillary  waves  present.  The  data  In  Table  6  clearly  Illustrate 
the  wide  range  of  c°  values  present  In  data  collected  over  Arctic 
regions. 


3 

OENSITY  DISTRIBUTIONS  ON  SAR  DATA 


The  distribution  of  the  recorded  Intensities  on  synthetic  aperture 
radar  (SAR)  Imagery  Is  the  result  of  radar  fading  (resulting  In  SAR 
Image  speckle)  and  the  distribution  of  the  scatterers  within  the  area 
Imaged  by  the  SAR.  The  radar  fading  dominates  the  spatial  patterns 
present  on  an  Image  unless  steps  are  taken  to  reduce  Its  Influence. 
Two  techniques  are  commonly  used  to  reduce  SAR  Image  speckle.  Non¬ 
coherent  Integration  or  multiple-looking  Is  achieved  by  processing 
separate  portions  of  the  SAR  bandwidth  Independently,  and  then  averag¬ 
ing  these  separate  sub-images  to  form  the  desired  output  Image 
(Porcello,  et  al.,  1976).  The  second  technique  simply  Involves  aver¬ 
aging  a  number  of  pixels  together  after  the  SAR  Image  has  been  formed. 
Both  of  these  methods  reduce  the  speckle  or  fading  In  a  SAR  Image  at 
the  expense  of  spatial  resolution.  These  methods  also  reduce  the 
amount  of  textural  Information  available  from  tne  SAR  data. 

To  better  understand  the  characteristics  of  the  Intensity  distri¬ 
butions  on  SAR  data  collected  over  distributed  targets,  an  analysis 
was  performed  where  the  observed  Intensity  distributions  from  airborne 
SAR  Imagery  were  compared  to  those  predicted  by  three  mathematical 
functions.  The  functions  evaluated  were:  the  Gamma  distribution,  the 
Inverse  Gaussian  distribution,  and  the  Log  Normal  distribution. 

The  SAR  data  from  each  test  site  was  first  digitally  processed 
Into  Imagery,  and  corrected  for  antenna  gain  and  range  fall -off  varia¬ 
tions.  No  multiple-looking  was  performed  during  the  processing  of  the 
SAR  data. 

An  80  by  80  pixel  subset  was  extracted  from  each  test  site  listed 
In  Table  7  for  each  SAR  channel.  Every  other  pixel  was  discarded  from 
this  Initial  sample  so  that  each  pixel  was  Independent  from  Its 
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TABLE  7 

SAR  IMAGE  INTENSITY  DISTRIBUTION  STUDY  TEST  SITES 


SAR 

TEST 

MISSION 

SITE 

U.S.G.S 

G1 

G2 

G3 

G4 

G5 

G6 

SARSEX 

SI 

S2 

S3 

S4 

S5 

S6 

S7 

MIZEX 

Ml 

M2 

DESCRIPTION 

Grass  Field 
3  Year  Old  Pine 
30  Year  Old  Pine 
50  year  Old  Pine 
Mixed  Oaks 

Sweetgum/Yellow  Popular 

Sea  Surface,  7  kt  wind 
Sea  Surface,  16  kt  wind 
Grass  Field 
Potato  Field 
Unharvested  Corn  Field 
Harvested  Corn  Field 
Mixed  Deciduous  Forest 

First  Year  Ice 
Second  Year  Ice 
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neighbors.  From  these  original  samples,  eight  additional  data  sets 
were  generated  by  spatially  averaging  between  2  and  9  adjacent  pixels. 
Thus,  for  each  original  80  by  80  pixel  data  set,  we  now  have  nine 
subsets,  with  the  number  of  Independent  samples,  N,  equal  to  1,  2,  3, 
4,  5,  6,  7,  8,  and  9. 

3.1  SAR  IMAGE  INTENSITY  MODELS 

The  Gamma  distribution  can  be  formulated  as 

f (x)  =  [(x/b)c-l  e-x/b]/[br(c)]  (3) 

where  x  Is  a  given  value  within  the  population,  r (c)  Is  the  Gamma  dis¬ 
tribution  evaluated  at  c,  and  b  and  c  are  defined  as 


b  *  s2/7 

(4) 

c  =  72/s2 

(5) 

where  7  Is  the  mean  value  for  all  x's  In  the  sample  and  s  Is  their 
standard  deviation.  The  Gamma  distribution  Is  used  to  model  radar 
fading  for  power  measurements  where  no  background  variation  Is  pre¬ 
sent.  In  this  case,  c  *  N,  where  N  Is  the  number  of  Independent  sam¬ 
ples  averaged  to  form  the  radar  data  set. 

The  Inverse  Gaussian  distribution  can  be  formulated  as 

(x)  “  {X/[2m(x-a)3]}l/2  •  exp{-  [x-<r-u]2/[2u2(x-a)]},  (6) 


where 


*  ■  *1  -  (x  -  xj)3/(2s2  log  n), 


(7) 
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xj  Is  the  minimum  value  of  x  In  the  sample, 

u  ■  J  -  a,  (8) 

and 

X  «  [(1/n)  E  (xi  -a)-1  -  (l/u)]-l.  (9) 

Eqs.  (7-9)  use  the  sample  data  set  to  generate  the  necessary  param¬ 
eters  for  the  Inverse  Gaussian  distribution  In  Eq.  (6). 

The  Log  Normal  distribution  can  be  formulated  as 

f(x)  ■  {l/[(x-a)a  sqrt(2f)]}exp{-0.5[ln((x-a)/m)]/a)2},  (10) 

where 

a  •  standard  deviation  (In  x),  (11) 

m  -  eu,  (12) 

with 

U  *  mean  (In  x).  (13) 

Again,  Eqs.  (11-13)  use  the  data  set  to  generate  the  necessary  param¬ 
eters  for  the  Log  Normal  Distribution  described  by  Eq.  (10). 

The  parameters  defined  by  Eqs.  (4,5,7-9,11-13)  are  maximum  likeli¬ 
hood  estimates  described  by  Cheng  and  Amin  (1981)  and  Folks  and 
Chhlkara  (1978). 

In  order  to  determine  whether  or  not  a  derived  Gamma,  Inverse 
Gaussian  or  Log  Normal  Distribution  matched  the  sample  population,  a 
non-parametrlc  statistical  test  was  applied.  The  type  selected  for 
this  analysis  was  the  Kolmogorov  test  described  by  Conover  (1981). 


This  test  operates  by  comparing  the  maximum  vertical  separation  be¬ 
tween  the  sample  distribution  and  a  candidate  distribution  [e.g.,  In 
our  case,  Eqs.  (2,5,9)]  with  a  Kolmogorov  test  statistic  (K).  The 
Kolmogorov  test  statistic  Is  based  on  the  number  of  points  In  the 
sample  population  (N)  and  the  level  of  significance.  For  a  95X  level 
of  significance,  this  statistic  Is 

K  =  1.36/N1/2  (14) 

This  test  statistic  Is  used  by  comparing  the  cumulative  distribu¬ 
tion  functions  of  the  candidate  and  sample  distributions.  Figure  16 
presents  a  plot  of  the  CDF  for  test  site  S3,  along  with  the  Inverse 
Gaussian  and  Log  Normal  distribution  fits  for  this  data  set. 
Figure  17  presents  plots  of  the  differences  between  the  actual  data 
and  these  two  distributions.  The  K-parameter  for  this  example  Is 
+  0.04.  Thus,  we  can  see  that  both  the  Inverse  Gaussian  and  Log 
Normal  distributions  adequately  model  the  observed  Intensity 
distributions  for  this  test  site. 

Cumulative  distribution  functions  for  all  three  speckle  models 
were  generated  for  all  test  sites  for  sample  numbers  N  «  1  to  9  using 
the  maximum  likelihood  estimates.  The  predicted  CDFs  were  then  com¬ 
pared  to  the  CDFs  from  the  actual  SAR  data. 

Surprisingly,  In  all  cases,  all  three  speckle  models  generated 
distributions  which  matched  the  SAR  data  using  the  Kolmogorov  goodness 
of  fit  test.  Overall,  the  best  fit  was  achieved  using  the  Gamma  dis¬ 
tribution.  In  the  next  section,  we  will  discuss  how  the  relationship 
between  N  and  the  c-parameter  generated  using  Eq.  (5)  can  be  used  to 
estimate  the  texture  within  the  SAR  scene. 

3.2  TEXTURE  MODEL 

The  variability  In  the  Intensities  recorded  on  a  SAR  Image  Is  due 
to  two  factors:  (1)  radar  fading  or  speckle;  and  (2)  the  natural 
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Figure  16.  Cumulative  Distribution  Function  (CDF)  for  Test 
Site  S-3  with  Inverse  Gaussian  and  Log  Normal 
Fits 
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DIFFERENCE 


2-MAR-3b 


KOLMOGOROFF  DIFF  ERENCE  FOTT  GRASS  1 

FOR  THE  INVERSE  GAIJSSIttN 


(a)  Inverse  Gaussian 


Figure  17.  Kolmogoroff  Difference  Plots  for  Test  Site  S3 
[Figures  (a)  and  ( b ) ] 
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variability  or  spatial  distribution  of  the  scatterers  within  the  scene 
resulting  In  the  radar  backscatter.  The  distribution  of  the  SAR  Image 
Intensity  can  be  accounted  for  by  using  a  multiplicative  model  which 
takes  Into  account  both  radar  fading  and  the  texture  present  In  the 
data.  This  variability  can  be  modeled  after  Ulaby,  et  al.  (1986)  as 

Plj  =  Pi  Ti(J)  FN(j)  (15) 


where 

Pm  Is  the  Image  Intensity  of  the  Jth  pixel  In  the  1th 
field, 

Pi  Is  the  mean  Image  Intensity  of  the  1th  field, 

Ti(J)  Is  the  texture  random  variable  accounting  for  the 
natural  variability  In  the  1th  field,  and 

Ffl(j)  *s  the  speckle  variable  accounting  for  signal  fading 

(this  fading  Is  characterized  using  a  r2  distribution 
with  2N  degrees  of  freedom). 

The  texture  random  variable,  T1(j),  In  Eq.  (15)  Is  dependent  on 
the  scattering  characteristics  of  the  surface  resulting  In  the  radar 
backscatter,  and  thus  should  be  considered  along  with  the  measured 
power  Intensity,  P-j j ,  when  developing  radar  scattering  models. 

Figure  18  presents  a  scatter  plot  of  the  Gamma  distribution  c- 
parameter  versus  the  number  of  Independent  samples,  N,  for  the  3  year 
loblolly  pine  field  (test  site  G3)  for  each  radar  frequency.  Also 
presented  In  Figure  18  are  the  least-squares  regression  lines  [12] 
which  best  fit  the  data.  Clearly,  a  strong  linear  relationship  exists 
between  the  c-parameter  and  N. 

The  variability  within  an  Image  due  to  the  texture  (st)  of  the 
target  resulting  In  the  radar  backscatter  can  be  expressed  after 
Ulaby,  et  al .  (1986)  as 
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St2  -  [(s/7)2  -  1/N]/[1  +  (1/N)]. 


(16) 


Combining  Eqs.  (16)  and  (5)  results  In 

st2  «  [(1/c)  -  (1/N) ] / [1  +  (1/N)].  (17) 

Note  that  In  order  for  Eq.  (17)  to  be  valid,  c  always  has  to  be 
less  than  or  equal  to  N.  Since  our  analyses  Indicate  that  c  Is  a 
linear  function  of  N,  l.e., 

c  -  AN  (18) 

where  A  Is  the  slope  of  the  regression  line  determined  by  least 
squares  techniques,  Eq.  (17)  can  be  expressed  as 

st2  *  [1  -A] /[AN  +  A],  (19) 

From  Eq.  (19),  several  observations  can  be  made.  First,  If  the 
slope  of  the  regression  line,  A,  Is  1.0,  then  st2  Is  0.  This  Is  con¬ 
sistent  with  our  earlier  observation  that  If  c  »  N,  then  the  Gamma 
distribution  describes  a  pure  speckle  background  with  no  scene  varia¬ 
bility.  Second,  for  any  N,  as  A  decreases,  st2  increases.  And  final¬ 
ly,  for  any  c,  as  N  Increases,  st2  decreases  (l.e.,  as  the  number  of 
pixels  averaged  Increases,  the  texture  within  the  scene  decreases). 
The  A-parameter  for  all  the  test  sites  examined  during  this  analysis 
are  summarized  In  Table  8. 
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TABLE  8 

SUMMARY  OF  A-PARAMETERS  FOR  STUDY  TEST  SITES 


TEST  SITE 

DESCRIPTION 

X-BAND 

C-BAND 

L-BAND 

61 

Grass  Field 

1.00 

0.80 

0.61 

G2 

3  Year  Old  Pine 

0.97 

0.83 

0.51 

G3 

30  Year  Old  Pine 

0.70 

0.77 

0.76 

G4 

50  Year  Old  Pine 

0.79 

0.72 

0.79 

G5 

Mixed  Oaks 

0.82 

0.67 

0.72 

G6 

Sweetgum/Yellow  Poplar 

0.92 

0.67 

0.72 

SI 

Ocean  Surface,  5  kt  wind 

0.83 

0.69 

S2 

Ocean  Surface,  15  kt  wind 

0.77 

0.87 

S3 

Grass  Field 

0.83 

0.61 

S4 

Potato  Field 

0.89 

0.C2 

S5 

Unharvested  Corn  Field 

0.59 

0.47 

S6 

Harvested  Corn  Field 

0.55 

0.92 

S7 

Mixed  Deciduous  Forest 

0.59 

0.64 

Ml 

First  Year  Ice 

0.70 

0.86 

M2 

Second  Year  Ic 

0.84 

0.73 

62 


4 

COMPARISON  OF  SCATTERING  COEFFICIENTS 
TO  SURFACE  ROUGHNESS  PARAMETERS 

Thus  far  In  this  report,  we  have  concentrated  on  presenting  means 
to  extract  L-band  clutter  statistics  from  SAR  Imagery.  These  tech¬ 
niques  have  Included  extraction  of  a0  values  from  SAR  data  as  well  as 
mathematically  modelling  the  distributions  of  the  o°  values.  Another 
goal  of  this  program  was  to  collect  Information  on  the  surface  rough¬ 
ness  characteristics  of  the  test  sites  being  examined,  and  to  corre¬ 
late  the  roughness  measurements  to  the  a0  values.  These  comparisons 
will  be  presented  In  this  chapter. 

In  Section  4.1,  we  present  scattering  models  which  utilize  surface 
roughness  parameters  to  estimate  radar  scattering  from  different  types 
of  fields.  In  this  analysis,  we  utilize  the  roughness  measurements 
collected  during  this  program  (summarized  In  Larson,  et  al.  1986).  In 
Section  4.2,  we  present  a  statistical  comparison  between  the  o°  from 
the  Ouke  Forest  Test  Sites  and  the  forest  parameters  for  the  22  test 
sites.  Finally,  In  Section  4.3,  we  discuss  the  Implications  of  the 
texture  measures  derived  from  the  Gamma  distribution  with  respect  to 
the  scattering  properties  of  the  scene. 

4.1  MODELLING  OF  RADAR  CROSS  SECTION  (ff°) 

Three  scattering  models  were  evaluated  for  the  describing  a0  val¬ 
ues  observed  from  several  of  the  Long  Island  Test  Sites.  These  models 
were  exercised  using  the  surface  roughness  measurements  collected 
colncldently  with  the  SAR  overflights.  Four  fields  were  used  in  this 
analysis:  Site  M  (smooth,  bare  field);  Site  H  (corn  stubble);  Site  K 
(grass  field);  and  Site  RF3  (grass  field). 

4.1.1  SCATTERING  MODELS 

The  scattering  models  utilized  In  the  present  study  were:  (1)  the 
Barrlck  and  Peake  (1967)  model  for  slightly  rough  surfaces;  (2)  the 


Barrlck  (1968)  model  for  slightly  rough  surfaces  based  and  Gaussian 
surface  height  statistics;  and  (3)  the  Rayleigh  scattering  model 
defined  by  Rice  (1951). 

The  two  scale  models  combine  the  backscatter  from  a  slightly  rough 
surface  with  the  scattering  predicted  by  specular  theory.  The  back¬ 
scatter  coefficient  as  derived  by  Barrlck  and  Peake  (BP)  for  a  slight¬ 
ly  rough  surface  Is  given  by 

ffxx  m  ^  ko4  cos^f j  l®xx^  W(-2k0s1n  J j)  (20) 


where 


oj  -  Incident  angle 

W(*i)  -  surface  height  spectrum 

a($l)  -  scattering  coefficient 

xx  -  transmitter/ receiver  polarizations 

Both  the  surface  height  spectrum  and  scattering  coefficient  have 
been  derived  from  ground  measurements  and  used  to  calculate  a  value 
for  a.  The  second  part  of  the  two  scale  model  uses  the  scattering 
coefficient  predicted  from  specular  theory  as  developed  by  Barrlck 
(1968).  The  scattering  coefficient,  derived  by  Barrlck  and  based  on 
Gaussian  surface  height  statistics,  Is  given  by: 

*vv  ■  *hh  *  [(sec4  0|)/$2]  exp  [(tan2  0i)/S2] IR(o) I2  (21) 

where  S2  ■  mean  square  of  surface  slope 

R(o)  -  normal,  0-0,  reflection  coefficient 

Values  for  the  parameters  In  the  BP  slightly  rough  model  are  derived 
from  surface  roughness  measurements  and  dielectric  constant 
measurements. 
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The  Rayleigh  scattering  model  as  extended  by  Rice  (Cosgrlff,  1960) 
Is  of  value  for  the  slightly  rough  surfaces  considered  In  this  Inves¬ 
tigation.  The  scattering  coefficient  o0  Is  given  In  Equation  22  for 
the  Rice  model . 

a0  *  8  cos4^  fi2kQ4J  *p(h)  J0(2kQh  sin  0^)  r  dr  (22) 

where  hz  ■  mean  square  surface  roughness,  and 

r)  *  surface  height  autocorrelation  function 

The  Barrlck  and  Peake  and  Rayleigh  (Rice)  scattering  models  were 
exercised  using  the  surface  roughness  and  moisture  parameters  listed 
In  Table  9  and  summarized  In  Figures  19  and  20.  Both  of  these  models 
were  exercised  for  VV  and  HH-polarlzatlons.  The  results  of  the  models 
are  presented  In  Figures  21  and  22.  Also  presented  In  the  plots  are 
the  SAR-observed  c°  values  for  the  HH-polarlzed  Imagery.  Finally, 
Figure  23  presents  a  comparison  of  the  HH-polarlzed  data. 

Overall,  we  can  see  that  the  Rayleigh  (Rice)  model  predicts  lower 
o°  values  than  the  Peake  and  Barrlck  model.  Both  models  predict  the 
observed  decrease  In  L-band  o°  as  a  function  of  Incidence  angle.  The 
Peake  and  Barrlck  model  results  match  the  observed  results  quite  well, 
with  the  model  overpredicting  the  SAR-observed  values  somewhat  for 
test  sites  K  and  RF3.  Contributions  to  the  scattered  field  from  the 
slightly  rough  component  of  the  BP  model  dominated  the  true  scale 
model  results.  The  only  test  site  where  the  Rayleigh  (Rice)  model 
comes  close  to  matching  the  SAR-observed  o°  values  Is  for  test  site  H. 
In  all  other  cases,  the  model  significantly  underpredicted  the 
observed  results. 
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surface  Profile  Field  M  (EW) 


DISTANCE  (IN) 

Figure  19.  Surface  Height  Distribution  Plots  for  Long 
Island  Test  Sites  [Figures  (a)  through  (d)) 
(a)  Site  M  -  Smooth,  Bare  Field 
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Figure  19  <b).  Site  H  -  Corn  Stubble 
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Figure  19  (d).  Site  RF3  -  Grass  Field 
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Correlation  of  Surface  Profile  Field  M  (NS)  correlation  of  Surface  Profile  Field  M  (EW) 
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Figure  21.  Results  of  Peake  and  Barrick  Model  Predicted 
versus  SAR-Derived  Values  of  a°  [Figures  (a) 
through  (d) ] 

(a)  Site  M  -  Smooth,  Bare  Field 
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Figure  21  (b).  Site  H  -  Corn  Stubble 
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Figure  21  (d).  Site  RF3  -  Grass  Field 
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Figure  22.  Results  of  Rayliegh  (Rice)  Model  Predicted 

versus  SAR-Derived  Values  of  O0  [Figures  (a) 
through  (d)] 

(a)  Site  M  -  Smooth,  Bare  Field 
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Figure  22  (b).  Site  H  -  Corn  Stubble 
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Figure  23.  Comparison  of  Modeled  versus  SAR-Observed  Values 
of  a° .  (Figures  (a)  through  (c)] 

(a)  SAR-Derived  Values 
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Figure  23  (b).  Peake  and  Barrick  Model 
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4.2  STATISTICAL  COMPARISONS 

The  next  step  In  the  analysis  was  to  correlate  the  radar  scatter¬ 
ing  coefficients  listed  In  Table  5  with  the  forest  canopy  parameters 
presented  In  Table  3.  The  test  sites  were  divided  Into  two  groups 
(hardwood  stands  and  conifer  stands)  for  these  correlations.  Simple 
linear  correlations  (Draper  and  Smith,  1966).  Indicated  a  significant 
relationship  (at  the  0.90  level  of  significance)  between  the  L-band 
scattering  coefficients  and  coniferous  tree  height  (  -  -  0.65)  and 
hardwood  stems  per  acre  (R  ■  0.69).  Multiple  linear  regression  anal¬ 
ysis  (Draper  and  Smith,  1966)  of  the  data  showed  a  significant  (at  the 
0.90  level)  correlation  between  the  C-band  scattering  coefficients  and 
deciduous  DBH,  basal  area  and  tree  height  (R  «0  .92)  and  the  L-band 
scattering  coefficients  and  coniferous  DBH  and  tree  height  (R  *  .90). 

These  statistical  correlations  should  be  beneficial  In  more  deter¬ 
ministic  modeling  efforts.  A  first  order  model  to  describe  the  forest 
radar  cross  section,  of,  can  be  expressed  after  Ulaby,  et  al .  (1982) 
as 

of  -  Tc2Ob  +  +  Tu2 (os  +  a$t)]  (23) 

where  Tc  Is  the  transmissivity  of  the  overstory  canopy  layer, 

Tu  Is  the  transmissivity  of  the  understory  canopy  layer, 

ob  Is  the  cross  section  of  the  branches  In  the  overstory, 

<*t  Is  the  cross  section  of  the  tree  trunks  In  the 

overstory  and  understory, 
cs  Is  the  cross  section  of  the  soil  layer,  and 
o$t  Is  the  soll/trunk  Interaction  scattering  coefficient. 

The  correlations  found  In  this  study  Indicate  that  at  C-  and  L- 
band,  a  significant  portion  of  the  radar  backscatter  Is  explained  by 
parameters  which  describe  the  dimensions  and  numbers  of  tree  trunks  In 


the  forest  stand.  These  results  Indicate  that  modeling  efforts  at  C- 
and  L-band  at  vertical  polarization  should  Initially  concentrate  on 
those  parts  of  scattering  model  described  In  Eq.  (23)  which  Involve 
the  o\  and  <7$t  terms. 

4.3  SCENE  TEXTURE  AS  A  FUNCTION  OF  SCENE  COMPOSITION 

Examining  the  scene  texture  as  described  by  the  A-parameter  pre¬ 
sented  in  Chapter  3  may  lend  Insight  into  the  processes  resulting  in 
the  radar  backscatter  from  the  different  test  sites.  For  this  study, 
six  test  sites  from  the  Duke  Forest  test  site  were  selected  (see 
Figure  13):  B9  (grass  field);  B8  (3  year  pine);  B19  (30  year  pine); 
BIO  (50  year  pine);  B7  (mixed  oak)  and  B5  (sweetgum/yel low  poplar). 

The  nature  of  the  scattering  will  affect  the  surface  texture  In 
two  ways.  First,  for  single-bounce,  surface  scattering,  the  spatial 
distribution  of  the  surface  scatterers  should  be  detected  directly  In 
the  observed  texture  pattern.  For  volume  scattering,  the  spatial 
patterns  of  the  scatters  will  not  be  as  clear  because  of  the  multiple- 
bounces  within  a  volume-scattering  medium  will  act  to  reduce  the  var¬ 
iation  of  the  pixel  Intensities  within  the  scene.  Forests  will  act  as 
a  complex  scene  with  respect  to  texture  because,  depending  on  wave¬ 
length,  parts  of  the  canopy  will  act  as  volume  scatterers  while  other 
parts  will  result  In  direct  radar  backscatter. 

Figure  24  presents  plots  of  the  regression  lines  of  the  c-  para¬ 
meter  as  a  function  of  N  for  the  six  test  sites.  These  plots  compare 
the  X- ,  C-  and  L-band  data.  Figure  24,  from  the  grass  field  (see 
Figure  13g),  contains  the  results  which  are  most  readily  Interpreted. 
At  X-band,  the  scattering  Is  most  likely  occurring  from  volume  scat¬ 
tering  from  the  grass  Itself.  Since  grass  Is  very  homogeneously  dis¬ 
tributed  throughout  the  scene  In  this  case,  with  no  bare  ground  being 
visible,  the  scattering  Is  almost  exponential  In  nature 
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Figure  24.  Regression  Lines  of  the  C-Parameter  as  a 

Function  of  N  (the  number  of  independent  pixels 
averaged) 
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(l.e.,  A  =  1.00).  The  opposite  extreme  from  the  volume  scattering  In 
the  X-band  case  Is  the  L-band  scattering,  which  Is  more  than  likely 
almost  entirely  from  the  soil  underneath  the  grass.  This  surface 
scattering  results  In  rougher  texture  In  the  L-band  Image.  At  C-band, 
we  would  expect  a  combination  of  volume  and  surface  scattering,  re¬ 
sulting  in  a  texture  somewhere  between  those  observed  at  X-  and  L- 
band.  This  Is  Indeed  the  case. 

The  next  most  complex  srattering  case  Is  represented  by  the  3  year 
pine  stand  (see  Figure  13a).  At  L-band,  the  branches  of  the  small 
pine  trees  are  probably  not  large  enough  to  influence  the  scattering 
of  the  vertical ly-polarizcd  waves,  but  the  trunks  of  these  trees  are 
probably  large  enough  to  result  In  additional  scattering  from  the 
scene.  This  additional  scattering  Is  Indicated  by  the  higher  L-band 
a0  value  for  the  3  year  pine  stand  relative  to  the  grass  field,  and 
also  results  In  more  texture  In  the  scene,  as  the  decrease  in  the  A 
term  In  Table  8  and  Figure  2 4  Indicates.  At  X-  and  C-bands,  there  Is 
almost  no  change  In  Image  texture  compared  to  the  grass  fields.  Indi¬ 
cating  that  the  vertically-polarized  electromagnetic  energy  at  these 
radar  frequencies  Is  not  being  greatly  Influenced  by  the  presence  of 
the  small  pine  trees  In  this  field. 

While  simple,  first-order  scattering  models  can  be  used  to  explain 
the  differences  In  observed  texture  for  the  grass  and  three  year  pine 
test  sites,  more  complex  models  are  needed  for  the  test  sites  which 
contain  a  fully-developed  forest  canopy.  The  only  mature  forest  test 
site  where  any  significant  variation  In  scene  texture  exists  Is  for 
the  sweetgum/yel low  poplar  stand.  Two  questions  arise  when  examining 
the  textures  for  the  mature  forest  stands.  First,  why  are  there  no 
variations  In  texture  as  a  function  of  radar  frequency  for  the  30  year 
pine,  50  year  pine  and  mixed  oak  stands?  And  second,  why  does  the 
Image  texture  at  X-  and  I. -bands  for  the  sweetgum/yel low  poplar  stands 
vary  from  the  other  stands?  The  answers  to  these  two  questions  depend 
to  a  great  degree  on  the  sources  of  the  scattering  at  the  different 
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radar  frequencies.  Radar  scattering  models  from  forest  canopies  have 
not  been  developed  to  the  point  where  these  questions  can  be  ad¬ 
dressed.  One  possible  explanation  lies  In  the  density  of  the  trees 
within  the  four  mature  forest  stands.  Prom  Table  5,  we  can  see  that 
the  sweetgum/yellow  poplar  stand  had  a  lower  tree  density  (470  stems/ 
ha)  than  any  of  the  other  three  sites  (620  to  963  stems/ha).  The 
lower  tree  density  may  result  In  a  discontinuous  canopy  profile  for 
this  stand  (versus  a  continuous  canopy  profile  for  the  others),  caus¬ 
ing  a  different  spatial  pattern  In  the  scene  scatterers,  which  Is 
detected  as  a  texture  difference  In  the  SAR  Imagery. 
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5 

CONCLUSIONS  ANO  RECOMMENDATIONS 


A  goal  of  research  sponsored  by  the  Rome  Air  Development  Center 
(RADC)  Is  to  develop  a  terrain  scattering  model  to  predict  terrain 
scattering  properties  at  L-band.  In  order  to  achieve  this  goal  there 
has  been  a  need  to: 

1.  Collect  a  data  base  of  radar  cross-section  signatures  (both 
the  mean  and  the  histogram)  for  various  terrain  types, 

2.  Collect  a  coincident  (with  [1]  above)  set  of  surface  scatter¬ 
ing  characteristics  (l.e.,  surface  roughness  and  dielectric 
property  measurements)  for  the  various  terrain  types,  and 

3.  Develop  a  model  (deterministic  and  statistical)  which  predicts 
the  radar  cross-section  signature  as  a  function  of  the  surface 
scattering  characteristics. 

Research  has  been  conducted  at  ERIM  for  RADC  over  the  past  two 
years  to  further  develop  L-band  terrain  scattering  models  (Larson, 
et  al.,  1986;  Kaslschke  1987).  Under  this  program,  a  set  of  airborne, 
L-band  SAR  data  collected  over  several  test  areas  has  been  analyzed. 
L-band  cross-section  signatures  of  various  terrain  types  have  been 
extracted  from  the  SAR  data  using  ERIM's  digital  SAR  Image  calibration 
algorithm  (see  Larson,  et  al.,  1986).  These  signatures  Include  both 
an  average  radar  cross-section  value  (o°)  as  well  as  a  radar  cross- 
section  histogram  for  each  terrain  type. 

A  selected  set  of  terrain  property  measurements  were  made  for  a 
limited  set  of  test  sites  during  each  SAR  overflight.  From  each  site 
the  following  Information  or  data  were  obtained:  (1)  soil  samples 
were  obtained,  from  which  dielectric  property,  soil  moisture  and  soil 
density  measures  were  obtained;  (2)  surface  photographs  were  taken  to 
document  the  general  state  of  the  vegetation;  (3)  the  height  and  den¬ 
sity  of  the  vegetation  were  noted,  and  measured  when  ever  possible; 
and  (4)  surface  height  profiles  were  obtained,  from  which  additional 
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surface  roughness  properties  were  derived,  Including:  (a)  spectrum  of 
the  surface  height,  (b)  slope  distribution,  (c)  radius  of  curvature, 
and  (d)  autocorrelation  of  the  surface  heights  (see  Larson,  et  al . , 
1985). 

Results  obtained  and  reported  under  the  present  research  effort 
Include  the  correlation  of  the  cross-section  signatures  with  surface 
parameter  measures  via  scattering  models.  These  models  are  either 
deterministic  or  statistical  In  nature. 

The  deterministic  models  are  used  to  correlate  the  average  radar 
cross-section  with  the  surface  scattering  properties.  The  second 
modeling  effort  has  utilized  methods  for  describing  the  statistical 
distribution  of  the  radar  cross-section  values  for  each  terrain  type. 
The  cumulative  distribution  model  Is  compared  to  the  actual  radar 
derived  distribution  via  non-parametrlc  statistical  tests. 

With  the  exception  of  Forest  sites,  all  of  the  test  fields  for 
which  surface  measurements  were  obtained  are  classified  as  slightly 
rough.  The  appropriate  surface  parameter  descriptor,  as  derived  from 
surface  height  and  dielectric  measurements,  are  utilized  In  scattering 
models  to  derive  backscatterlng  coefficients.  Results  obtained  using 
the  Barrlck  and  Peake  model  for  slightly  rough  surfaces  provided  best 
agreement  with  o°  values  derived  from  calibrated  SAR  data  of  the  test 
fields.  Calculated  values  using  the  Rice  model  were  6  to  10  dB  larger 
than  measured  cross  section  values.  Histograms  of  o°  for  various  test 
fields  were  constructed  from  SAR  data.  An  Inverse  Gaussian  distribu¬ 
tion  provided  the  best  fist  to  the  histograms  of  the  test  fields.  The 
Gamma  distribution  provide  best  fit  to  the  Forest  test  sites. 


The  limited  number  of  fields  and  range  of  roughness  values  pre¬ 
cluded  extensive  analysis;  however,  the  results  obtained  did  Indicate 
a  trend  of  Increasing  a0  with  Increasing  rms  surface  height  roughness 
scale  for  bare  fields.  Correlation  lengths  of  surface  roughness  did 
not  show  any  particular  correlation  with  o°. 
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Derived  values  of  aQ  have  been  compared  to  other  published  values 
as  obtained  from  the  University  of  Massachusetts  data  bank.  Using 
published  results  obtained  with  conditions  that  best  fit  our  derived 
values,  general  agreement  (within  2  to  5  dB)  were  noted.  Data  uti¬ 
lized  In  a  comparison  Included: 

•  A.J.  Blanchard,  et  al.,  IEEE  Geo.  Rem.  Sen.  GE20(1),  pp.  36- 
41,  1982;  and 

•  F.T.  Ulaby,  et  al.,  IEEE  Geo.  Rem.  Sen.  GE17(2),  pp.  33-40, 
1979. 

Most  data  were  obtained  with  relatively  large  values  of  Incident  angle 
so  that  comparisons  are  difficult  as  little  low  angle  results  are 
available. 

Results  obtained  have  shown  that  terrain  scattering  coefficients 
can  be  derived  from  calibrated  L-band  SAR.  In  addition,  the  large 
areal  coverage  capability  of  the  SAR  Is  used  to  obtain  data  from  which 
scattering  coefficients  from  a  large  number  of  sites  and  terrain  types 
are  derived.  Although  only  a  small  fraction  of  the  calibrated  L-band 
SAR  data  available  were  utilized  In  this  research  due  to  funding  limi¬ 
tations,  scattering  coefficients  have  been  obtained  from  a  wide  vari¬ 
ety  of  terrain  types  and  forests  areas.  Data  may  be  obtalred  for  a 
wider  range  of  Incident  angles,  being  limited  only  by  the  number  of 
SAR  data  collection  passes  made. 

Research  reported  In  this  report  has  focused  on  the  reduction  and 
analysis  of  L-band  data  from  several  programs  conducted  for  other  U.S. 
Government  sponsors.  Although  much  progress  was  made,  there  are  still 
fundamental  limitations  In  this  data  set.  These  limitations  are  with¬ 
in  two  areas:  (l)  the  number  of  test  sites  and  their  location  limited 
the  variation  In  surface  scattering  parameters  which  Is  needed  to 
fully  develop  L-band  terrain  scattering  models;  and  (2)  techniques 
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have  to  be  further  developed  to  measure  the  surface  scattering  proper¬ 
ties  of  terrains,  especially  tree  canopies.  The  following  are  specif¬ 
ic  recommendations  for  future  work  to  continue  and  address  the  above 
stated  goal  of  RADC. 

1.  Collect  L-band  radar  cross-section  coefficients  over  a  wide 
variety  of  surface-terrains  under  a  variety  of  conditions 
using  a  surface-based  scatterometer, 

2.  Utilize  available  L-band  SAR  data  and  where  necessary  collect 
additional  radar  cross-section  data  over  different  terrain 
types  as  required,  obtain  surface  measurements  Incident  with 
RCS  measurements, 

3.  Evaluate  new  techniques  to  measure  the  surface  characteristics 
of  various  terrain  types.  Including  the  use  of  a  laser  profil¬ 
ing  system  and  photographic  techniques,  and 

4.  Correlate  the  surface  scattering  measurements  to  the  L-band 
radar  scattering  data  via  deterministic  and/or  statistical 
scattering  models. 

It  Is  anticipated  that  arrangements  may  be  made  In  the  future  to 
collect  X-L-C-band  SAR  data  using  the  ERIM  SAR  which  Is  being  deployed 
on  a  U.S  Navy  P-3.  With  sponsor  approval,  the  data  from  this  system 
may  also  be  processed  to  generate  X-band,  L-band  and  C-band  terrain 
radar  cross-section  values.  It  Is  recommended  that  surface-based 
scatterometers  be  used  to  collect  data  at  sites  within  each  test  area 
coincident  with  SAR  data  collection.  The  SAR  and  scatterometer  data 
collections  should  be  coordinated  so  that  each  Instrument  Is  collect¬ 
ing  data  on  the  same  day.  The  scatterometer  can  be  configured  to 
collect  data  at  a  range  of  Incidence  angles. 

It  Is  further  recommended  that  at  each  test  site  during  SAR  and 
scatterometer  data  collection  period,  surface  scattering  data  be  col¬ 
lected,  Including  soil  samples,  surface  profiles  and  surface  photo¬ 
graphs.  One  site  should  be  established  as  an  intensive  test  site. 
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Within  this  site,  extensive  scatterometer  and  surface  scattering  mea¬ 
surements  should  be  obtained,  Including  surface  roughness  profiles 
utilizing  a  laser-profiling  device  such  as  what  has  been  developed  at 
ERIM. 

Verification  of  methods  for  the  characterization  of  surface  rough¬ 
ness.  The  use  of  recently  developed  laser  profilers  that  allow  large 
areas  to  be  profiled  with  accuracies  of  better  than  1  cm  should  be 
utilized  In  future  research.  Surface  contours  obtained  using  this 
system  can  be  correlated  with  samples  of  surface  roughness  obtained 
using  conventional  techniques  (Larson,  1986)  to  determine  when  surface 
roughness  profile  has  been  adequately  sampled  for  use  In  scattering 
models.  Since  the  measured  surface  profile  Is  used  to  derive  other 
surface  descriptors  for  use  In  scattering  models,  the  answer  to  this 
speculation  Is  critical  to  the  development  of  reliable  models. 

It  Is  recommended  that  the  North  Carolina  test  be  used  to 
further  develop  techniques  to  quantify  the  structural  characteristics 
of  the  forest  canopy,  and  to  relate  these  characteristics  to  the  radar 
scattering  signatures.  This  Is  recommended  because  of  the  extensive 
ground  truth  data  available  describing  the  entire  forest  test  site. 
The  surface  profiling  techniques  developed  for  low-lying  vegetation 
canopies  cannot  be  applied  to  forest  canopies,  and  different  ap¬ 
proaches  must  be  developed.  We  recommend  photogrammetrlc  techniques 
to  measure  the  canopy  closure  and  the  horizontal  and  vertical  surface 
area  of  the  different  forest  stands  within  the  experimental  forest 
west  of  Durham,  North  Carolina  (l.e.,  the  Duke  Forest).  These  mea¬ 
surements,  along  with  the  other  canopy  characteristics  which  are  on 
record  for  this  forest  (l.e.,  canopy  height,  stand  or  stem  density, 
age,  species  diversity,  and  basal  area  or  average  stem  diameter)  may 
be  correlated  with  the  L-band  radar  cross-section  data  collected  over 
this  area.  A  data  set  already  exists  where  SAR  data  were  collected 
when  there  were  no  leaves  on  the  trees.  It  Is  recommended  that  an 
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other  data  set  be  collected  when  leaves  are  on  the  trees  along  with 
coincident,  photographies  data  when  the  SAR  data  is  collected. 

Analysis  of  this  new  data  would  provide  the  radar  cross-section 
data  to  be  correlated  to  the  surface  scattering  data  using  statistical 
or  deterministic  scattering  models,  much  In  the  same  manner  as  was 
used  In  the  present  program.  These  results  may  then  be  used  to  add  to 
the  available  data  base  for  an  L-band  scattering  model,  and  to  pro¬ 
vide  additional  empirical  data  correlation  for  verification  of  the  L- 
band  scattering  model. 
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